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Abstract 

The lattice Boussinesq equation (BSQ) is a three-component difference-difference equa- 
tion defined on an elementary square of the 2D lattice, having 3D consistency. We write the 
equations in the Hirota bilinear form and construct their multisoliton solutions in terms of 
Casoratians, following the methodology in our previous papers. In the construction it turns 
out that instead of the usual discretization of the exponential as [{a + k)/{a — /c)]" we need 
two different terms [{a — ujk)/{a — fc)]" and [(a — ui'^k)/{a — fc)]", where a; is a cubic root of 
unity ^ 1. 

1 Introduction 

Among the integrable 2D difference equations one important class consists of those equations 
that are defined on an elementary quadrilateral and are multidimensional consistent. Here 
multidimensional consistency means that one can add a third dimension and extend the definition 
in a natural way from a quadrilateral to a cube, and that on this cube the maps are consistent. 

Most of the maps in this class are included in the ABS classification [1], (which is complete 
within the assumptions of symmetry and the tetrahedron property) . These lattice maps have one 
component for each lattice site, but there are also 3D-consistent multicomponent maps related 
to the Boussinesq (BSQ) equation [HEIHIIS]. 

In this paper we will derive multisoliton solutions to the lattice BSQ equation defined on the 
elementary square by the equations [5] 

(1.1a) 
(1.1b) 

0, (1.1c) 



Bl 


= w 


— UU + V = 


0, 


B2 


= w 


— UU + V = 


0, 


B3 


= w 


— uu + v — 


p-q 



u — u 



where we have used the standard shorthand notation, e.g., u = Un+i,rm v = Vn,m+i, and where 
p and q are parameters in the n and m directions, respectively. From equations (jl.lap . ()l.lbp 
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one finds, after taking the hat and tilde shifts, respectively 



u 



w 



V — V 

u — u^ 
uv — uv 
u — u 



(1.2a) 
(1.2b) 



Equations ()l.ip . (|1.2p are consistent on the elementary square. 

Three-dimensional consistency means that we can add a third direction, with parameter r. 
The new equations will be 



w — UU + V 
W — UU + V 



w — UU + V =0, 
r — 



u — u 
r — p 



=0, 



u — u 



(1.3a) 
(1.3b) 

(1.3c) 



where we have denoted the shift in the third direction by a bar. To this we should add the 
bar-tilde and bar-hat versions of (jl.2p . which can be derived from (jl.3ap . 

Consider now the cube of Figure [H where F stands for the three components {u,v,w). The 



F. 



F 



F 



F , 



F 



F 



Figure 1: The multi-dimensional consistency cube. 

initial values are given at F, F, F, F, but due to the linear equations (jl.lap . p".lbp . (jl.3ap we only 
need to give {u,v,w), {u,v), {u,v), {u,v). Given these values we can compute (u,v,w), {u,v,w) 
{u,v,W) using (|1.2p and its tilde-bar and hat-bar versions, along with (jl.lcp . (fr.3bp . (jl.3cp . After 

this there are three different ways to compute the remaining (n, v, w) but they all give the same 
values and the system is therefore SD-consistent. 

In the following construction we make active use the 3D-view with different interpretations 
for the bar shift. In Section [2] we construct first the background solution, then the one-soliton 
solution, and then in Section [3] we derive the Hirota bilinear form and propose and prove the 
formulae for A^-soliton solutions. The procedure is in principle similar to the one used in [6l[7l[8], 
but the multicomponent nature induces some new features. 

2 The background and one-soliton solutions 

The fixed-point idea proposed in [6j for the construction of the background solution (OSS) is 
that the background or "seed" solution should be a fixed point with respect to a shift in the 
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third direction, i.e., u = u,v = v,w = w. The relevant equations (the equations on the sides of 
the cube) are then obtained from (11.31) as 



W = V? — V, 



p—r 

u—u 
^ ^ I r—n 

w = uu — V + 



w = uu-v + ^, (2.1) 



where r is the parameter in the bar-direction and now plays the role of a parameter of the 
background solution. The equations are easy to solve and we find the OSS 



uq = an + bm + ci, (2.2a) 
VQ = \ul + \{a'^n + b^m + C2)+cz, (2.2b) 
wq = \ul- \{o?n + h^m + C2) - cz, (2.2c) 

where a, h are related to p, q by 

a^ = r-p, b^ = r-q, (2.3) 

and ci,C2,C3 are arbitrary constants. 

The one-soliton solution (ISS) is constructed with the same idea of using the 3D cube, but 
now u,v,w correspond to the background solution and u,v,W to the ISS. In more detail, we 
have the side equations 

w = uu — V, (2.4a) 

together with 

33 I ^ 3 3 (2 4b) 

nu-w + ^, \v = uu-w + '-^. 

These act as a Backlund transformation (BT) with k as the BT parameter related to the bar 
direction. 

In order to solve equations (j2.4p we take {u,v,w) = {uo,vo,wo) where the background 
solution is defined in (12.21). The ISS is written in the form 



{u,v,w) = {uo + x,vo + y,wo + z), (2.5) 

where {uq,vo,Wo) is the bar-shifted {uq,vo,wo), i.e., 

Uq = an + bm + k + ci, (2.6a) 

vo = ^ul + \{c?n b^m A; V C2) + C3, (2.6b) 

wq = ^Uq — ^{a^n + b^m + k"^ + C2) — C3. (2.6c) 



With these definitions we find from ()2.4ap that z = uqx. Thus we only need to solve for x,y, 
for which we have from (|2.4bp 



-up x+y ( y = -{vo+wo)x+uoy 

x—a+k ' / _ x—a+k ' (o 7) 

-up x+y I _ -{vo+wo)x+uoy ^ ' ' 

x-b+k ' L y x-b+k ' 
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This system can be linearized by taking {x,y) = (y, the result is 

$ = AT^, $ = M^-, (2.8a) 

where 

io -1 \ / % -1 \ 

^i>= \ H\ , N = \ '^o + wo -uo \ , M = Ivo + wo -uo . (2.8b) 

-1 a-kj \ -1 b-kj 

The matrices A^, M satisfy the integrability condition NM = MN. 

In order to construct the solutions it is useful to note that if we define 

(Uo{n,m)—ujk —1 0\ 

uo{n,m) -uj'^k -1 , (2.9) 

(-no(n,m) + A:)/(3A:2) l/{2,k^) l) 

then 

^a-ojk \ 

N = Q{n + l,m)-^ D{a)Q{n,m), where £'(a) = ( a-uj'^k , (2.10) 

a-k) 

and similarly M = Q{n,m + 1)~^ D(6) (5(^^, w-). Here a; is a cubic root of unity, cu ^ 1, i.e, 
+ cj + 1 = 0. Using these we can straightforwardly construct the solution as 

^(n, m) = Q{n, m)-^ D{a)''D{b)'^Q{0, 0)^(0, 0), (2.11) 

from which we find 

G = k{uj-l)[p'i{a-ujk)''{b-ujkr -Lo^p'^^ia-Lu^krib-Lo^kr], (2.12a) 
H = uog + k\u-l)[-uj^p°^{a-ujk)''{b-ujkr + p^2{a-uj^k)'^{b-uj^kr],{2.12h) 
F = pl{a-kY{b-kr + plia-ujkTib-uk)'^ + pl{a-u'^kT{b-0J^k)'^, (2.12c) 

where we have introduced new constants p^ in place of Gqo, -ffoo; -^oo- Using ()2.5p . ()2.6p we can 
recover the ISS as 

u''' = u, + k^-±^^^P^±^. (2.13a) 

1 + pi + P2 

v'^' = .0 + no A;l±i^^^l±i^ + k^ ^+^'Pi+^P2 ^ (2.13b) 

1 + Pl+ P2 1 + Pi + P2 

15-5 , , 1 + wpi + ^^2 iQ N 

w = Wo + tio K • (2.13c) 

1 + Pl + P2 



Here uq^vq^wq were defined in (12. 2p and 



^ (a -A;)" (6-A;)"^ pg ^ ^ 



where A; is the soliton parameter. 
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3 Bilinearization, Casoratians and A^-soliton solutions 
3.1 The main result 

The ISS (j2.13p is not quite sufficient for guessing the general structure for NSS, but after 
considering also 2SS (in the case = 0) we arrived to a solution in terms of Casoratians, 
constructed as follows: Given the multi-indexed function 

3 

il}j{n,m,l) = ^Qfl{5 - uj'kj)\a - uj'kjYih - uj'kj)"^ , (3.1) 

s=l 

we define the column vector 

ip{n,m,l) = {■ipi{n,in,l),tlj2in,m,l),- ■ ■ ,tp]y{n,mj))'^ , (3.2) 

and then the generic N x N Casorati matrix by combining columns with different shifts /j. The 
generic Casoratian is then the determinant 

Cn,m{i';{li}) = \'ip{n,m,li),'ip{n,m,l2), - ■ ■ ,ip{n,m,lN)\- (3.3) 

To describe such Casoratians we use the shorthand notation [9] in which only the shifts are 
given. Furthermore for consecutive sequences we use M = 0,1, . . . , M (this cannot be confused 
with the use of hat for shifts). 

Proposition 1. Multisoliton solutions to (11.11) are given by 

9 9 h 9 s 

U = Uq--, v = Vq-Uqj + -, W = Wq-Uq- + -, (3.4) 

where uq,vq,wq are given in ()2.2p and the functions f, g, h, s are given in terms of Casoratian^ 
f = \N^\, g = \N~^2,N\, h = \N^,N + 1\, s = \N^3,N -1,N\, (3.5) 
composed of ip given in (|3.1|) with 5 = 0. 

Here the size on the matrix indicates the number of solitons and the set {/cj}^^ provides 
the "velocity" parameters of the solitons, while the parameters gj^l are related to the locations 
of the solitons (by gauge invariance only their ratio is significant). 

In order to prove the above Proposition, we note that using ()3.4p as a dependent variable 
transformation we can bilinearize (|l.ip as 

Bi = J{h + ag) - g{g + af) + fs = 0, (3.6a) 
132 = f{h + bg) - g{g + bf) + fs = 0, (3.6b) 

B3 = fg-fg-{a- b){ff- fj) = 0, (3.6c) 

Ba = {a" + ab + b^)if7- If) + (a + 6)(/f - Jg) +}s + fh-gl = 0. (3.6d) 



If Af = 1 then / = \0\,g= \l\,h= \2\,s^ 0; and if iV = 2 then / = \0,l\,g= \0,2\,h= 10,3l,s = \1,2\ 
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In fact, the lattice BSQ equation (jl.ip can be written as 



m = Sl B2=Sl (u - 5)53 = "^"^ - '° - '^""^ + + "''""^ . (3,7) 
// // //// 

Thus we need to prove that the Casoratians ()3.5p solve the bilinear equations ()3.6p . This is given 
in Appendix |A] for the generic values of 6, although only 5 = is used for the lattice BSQ. The 
Casoratian proof suggested the 5^0 generalization along with some others, they are discussed 
next. 



3.2 The role of 5 

The solution to the lattice BSQ are obtained with (5 = in the matrix entry (j3.ip but the 
generalization 5 7^ is natural and we may ask about its meaning. It is important to note that 
the proof given in the Appendix can be carried out using only the following assumptions on the 
column vectors 

{a-5)ijj = V-?, (3.8a) 
(6-5)^ = V-?, (3.8b) 

7^ = V'-3# + 3(^V, (3.8c) 

where 7 is a diagonal matrix and the bar-shift of ip is defined by •0(n,m, /) = ip{n,m,l + 1). 
From these assumptions one can derive the bilinear equations 

Bi = f[h + {a- 5)g] - g[g + {a - 6)f] + fJ = 0, (3.9a) 

Bi = f[h + (6 - S)g] - g[g + (b - 5)f] +fs = 0, (3.9b) 

Bl = fg-fg+ (a - b){ff- //) = 0, (3.9c) 

Bi = {a^ + ab + b^)iff- Jf) + {a + b + 5){Jg - /f) + Js + fh-g1 = Q. (3.9d) 



We can now reverse the dependent variable transformation (13. 4p and construct from (13. 9p a 
generalized lattice BSQ equation 

w = uu — V + 5(u — u — a), (3.10a) 
w = uu — V + 6{u — u — b), (3.10b) 

w = uu — V -\ — — — 6(u — u — a — b). (3.10c) 

u — u 

Obviously, if we take 5 = 0, the above equations reduce to those of lattice BSQ. 
Now considering the Casoratian forms of f,g, h, s one can easily show that 

/(^) = /(O), 

g{6) = giO)+N5fiO), 

h{6) = h{0) + {N + l)6g{0)+N{N + l)/26^ f{0), 

s{S) = s{0) + {N -l)6g{0)+N{N -l)/26^ f{0), 
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where is the dimension of the matrix. From this result, which was derived from a particular 
form of the solution, we are led to the following: If we denote the (5-dependent functions in ()3.10|) 
by u',v',w' then the transformation 

u' = u, v' = V — S {u — uq), w' = w + 6 {u — uq), (3-11) 

converts ()3.10p into (jl.ip . Thus the effect of introducing 5 in (j3.ip can be undone by the "gauge" 
transformation (j3.1ip . 



3.3 Toeplitz generalization 

Another generalization in the construction of the solution is obtained if we replace (j3.8cp with 



(3.12) 



where T is some N x N matrix. Under this generalization proof of the bilinear equations 
proceeds as before, except for some details described in Appendix I A. 3[ 

First note that since the solutions are given in terms of determinants any matrix similar to 
F yields same solution as F. Thus it is sufficient to consider only different canonical forms of F. 

If F is a diagonal matrix 

F = Diag(7i,72,--- ,7Ar) (3.13a) 

where 

= 5^ - /c3, (3.13b) 

with distinct {kj}, then the entries in the Casoratian can be as in (j3.ip . Since = lu* (where 
* stands for complex conjugate), the condition for a real solution (coming from real ipj) is 



.(0) 



.(0)* 



with kj , gfl E 



Next suppose that F is a lower triangular matrix defined as 



^Nikl) = {7s,likl))NxN, ls,iiki) 



(s-l) 





s < I, 



(3.14) 



(3.15) 



where 71 is defined by ()3.13bp . In this case, the generic entry vector ip can be taken as 



ij = ^AsQsiki), 



s=l 



where 



Qsiki) = (QsAh), Qs,i{ki), ■■■ , Qs,N-i{ki)y 

Qs,j{ki 



Iff) 



gfliS - Lo'ki)\a - u)'kiY{h - uo'kiY 



and {As} are arbitrary A^th-order lower triangular Toeplitz matrices defined by 



a^j G C. 



NxN 





/ 
















\ 









































v 


as,N-i 


0'S,N-2 


as,N-2, ■ 


• as,i 


O-sfl 


/ 



(3.16a) 



(3.16b) 
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To get a real solution, in addition to the condition (j3.14p we also need A2 = and ^3 being 
real. Reduction to the solution of the lattice BSQ (jl.ip are obtained by taking 5 = in ()3.16bp . 

Note that in this case there is only one ki, so solutions corresponding to such a rjv(^i) are 
kind of limit solutions. In fact, let us start from the Casoratians with ipj defined in (|3.ip and set 
g^^l = gf'l for j = 2, 3, • • • ,N. First we replace the Casoratians /, g, h, s in the bilinear equations 

(13:9]) by f/K,g/K,h/K,s/K with K = ^'''[j-iy'' ■ Then, using L'Hospital rule we can 

take the limits kj ki from j = 2 to j = N step by step, and finally reach the bilinear equations 
()3.9p solved by Casoratians /, g, h, s with entry vector ip = Yll=i Qs(^i) where Qs{ki) is defined 
as (|3.16bp . The Toeplitz matrix As in p.l6ap can be obtained by a suitable redefinition of g^^l 

as a function sufficiently differentiable w.r.t ki and g^^l = gf'likj) for j = 2, 3, • • • , A^. 

From the above solution we can also derive rational solutions, by taking a particular choice 
of parameters in the limit ki — > 0. 

4 Conclusions 

In this paper we have bilinearized the lattice Boussinesq equation (II. ip and constructed its multi- 
soliton solutions in terms of Casoratians. The method, which relies heavily on 3D-consistency, 
is similar to the one used in [71 Ej , with the main difference that the analogue of the exponential 
factor now contains cubic roots of unity, see (j2.14p , allowing two different p-terms in the bilinear 
construction. Cubic roots of unity also enter in the continuum case, because the continuum 
Boussinesq equation is obtained as a three-reduction from the KP equation. 

After this work was completed we were informed about reference |12] where the authors 
bilinearized the lattice BSQ in its 9-point 1-component form using singularity confinement and 
then constructed its solutions in Casoratian form. Their final result is similar to the present 
except that it only contained two terms in the Casoratian entries (cf. (13. ip ). 
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A Proof of Casoratian solutions 

Usually in Casoratian proof the bilinear equation is reduced to a Laplace expansion of a 2A'^ x 2A'^ 
determinant that can be seen to be identically zero. The expansion is described as: 

Lemma 1. Suppose that B is an N x (N — 2) matrix and a, b, c, d are N 'th-order column 
vectors, then 

|B,a,b||B,c,d| - |B, a, c||B, b, d[ + [B, a, dI|B, b, c[ = 0. (A.l) 

This Lemma is also used to generate Casoratian equalities by means of which one can simplify 
Casoratian proofs. 

A.l Formulae for Casoratians 

In order to use the above Lemma we need various formulae for the shifts of the Casoratians 
f,g,h,s of (j3.5p with entries (j3.ip . as given below. These formulae can be derived using 
()3.8ap . ()3.8bp in the same way as in For convenience we introduce an up-shift operator 
E" hy ~ ^ _ 

Down shifts are denoted by E^, v = 1,2,3, respectively. 
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The basic shift formulae are 



[a. 



N~2 : 



[a. 



N-2 i 



\N-2 ; 



{a-b){a-6) 



N-2, 



\N-2 . 



(a - 6)(a - bY^-\h - b)''-' g + {a + b- 25) f 



jiV - 2,E^7p{N - 2)1, 




(A. 2a) 


\N^3,N-l,E^iP{N- 


2)1, 


(A.2b) 


\N^3,N,E^ij{N -2)1 




(A.2c) 


\N^3,ij{N -2),i;{N - 


-2)1, 


(A.2d) 


-2,i;{N-2),i^{N-2)\, 




(A.2e) 



(a -b){a- 5)^-^{b - 5^"-'' s + {a + b- 26)g + [(a - 6)'' + {a - 6){b - 6) + {b - 6)'']f 



\N-2 



iV y I\ \ 

^|ai,--- ,a.j-i, haj, aj+i,--- ,aAr| = f j[ai, 
.7 = 1 S = i ^ 



(A.3a) 



= -{b-5){a- Sy^-'f + {a-6){b- 6)'^^'/ + \N - 5,N - 3,N - 2, ^{N - 2),i^{N - 2)\. 

(A.2f) 

where /U = 1, 2, and ai = a, 02 = b. 

In order to apply Lemma [T] in the proof of the main result we need some further equalities 
(also derived using Lemma [TJ: 

f\N^5, N -3,N- 2, iPiN - 2),ij{N - 2) \ 

= -{b- 6)^-^f\N^5, N -3,N -2,N -1, i){N - 2)| 
+ {a- 5)^-'^f\'N^,N -3,N -2,N -l,i){N -2)1 
fllf^A, N - 2,7p{N - 2),7p{N - 2) | 

= - (6 - (5)^"V|^V^, - 2, iV - l,Tp{N - 2)1 
+ (a - 6)^-^f\lf^4, N-2,N-1, i!{N - 2)\ 
={b - 5)^- V [{a - 8)''~'s + (a - 5)^-^5 + (a - 5)^/ - / 
- (a - <5)^-2/ [{b - 5)^-2. + (6 - 5)^-1^ + (6 - 5)^/ - / 

g|iV^, A^ - 2,'4,{N - 2),^{N - 2)\ 

= -{b- 5)^-2^|]V^, AT - 2, AT, i;{N - 2)\ 
+ (a - 5)^-2/|]V^, AT - 2, AT, V'(Ar - 2)|, 
(a - 6)/g =/[/! + (a - 5)^] - /[/i + (6 - 5)gl 
(a - b){a - 6)^-\b - 5)^-^ fh = {a - 6)^'^ f\N^3, N + l,f{N - 2)\ 

- (6 - (5)^~2_^|]V^, Af + 1, ^(A^ - 2)|. 
For the next identity we also need the following Lemma: 
Lemma 2. 

,a7v|, (A.4) 



(A.3b) 

(A.3c) 
(A.3d) 

(A.3e) 
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where aj = (aij,-- - ,a]\ij)'^ and b = (6i,--- ,b]\f)'^ are N 'th-order column vectors, and haj 
stands for {biaij, • • • , bNONj)'^ ■ 

Then, noting that 

(<5 - ujkjf - 36{6 - ujkjf + 36^6 - ujkj) = {6- kjf - 36{6 - kjf + 36^{6 - kj) = jj, 
we get, using Lemma [2] and the following identity 



N N 

[(E7.)/]/=[(E7.)/]/ 



(A.5) 



j=i j=i 

the explicit result 

(a - 6)^'^f [|iV^, N -3,N -2,N - 1, -^(iV - 2) [ 

- |iV^, - 2, iV, 7p{N - 2) I + |]V^, + 1, ij{N - 2) 
+ {b-6f+'f + g-ib-5)f 

- 36{b - 6f-\s -h)- 36\b - 5f-^g 

- (6 - - 3, - 2, - 1, i^{N - 2) | 

- |]V^, A^ - 2, A^, i){N - 2) I + jiV^, N + 1, ^{N - 2) | 
+ (a-<5)^+V + 5-(a-<5)/ 

- 3^(a - 5)^-2 (s -h)- 36^ {a - Sf-^g 



(A.6) 



A. 2 Proof for bilinear equations (13 .Qp 

With these Casoratian formulae derived in Appendix lA.il we can prove bilinear equations ()3.9p . 



Proof for (|3.9ap : We consider the down-tilde-shifted B\, i.e., 

f[h + (a - 6)g\ - g[g + (a - 5)f] + fs = 0. 

Using (|A?2cl) . fOb]) and (jX2al) we have 

-(a - 5)''-^[f[h + (a - <5)5] - ffb + (a " -5)/] + Z^] 
= \lf^l\\lf^3,N,Tp{N - 2)1 - |]V^, A^lliV^, Af - 1,V'(A^ - 2)| 
+ |]V^,V'(Af-2)||]V^,Af- l,Af|, 

which is zero in the light of Lemma [U by taking 

B = (]V^), a = V(A^-2), h = ip{N-l), c = Tp{N), d = 7p{N-2). 

()3.9bp can be proved similarly. 



(A.7) 
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Proof for (j3.9cp : We prove it in its down-tilde-hat-shifted version: 

f[g + (a - 5)f] - f[g + (b - 5)f] - (a - b)ff = 0. (A.8) 

This can be verified by using first ()A.2ap . ()A.2bp . ()A.2dp and then Lemma [1] with B = {N — 3), 
a = tp{N - 2), b = ip{N - 1), c = ip{N - 2) and d = ^{N - 2). 

Proof for ()3.9dp : We first rewrite ()3.9dp in the following form 

f[s + {a + b + 6)g+ia^ + ab + b'^)f] - g\g + (a + b + 6)f] - (a^ + ab + b"^)/ f + fh = 0. (A.9) 

Then using (jA.2eP and (|A.2f|) we have 

(a - b){a - 6)^-\b - 6 f-^ x l.h.s.^EM 
= f\N^5,N-3,N-2,4:{N-2),i;{N -2)\ 
+ (35/ - g)\N^A, N -2,^Ij{N -2),i;{N -2)\ 
+ {a - b){a - 5)^-2(6 - 5)^-^(3^2^ _ ^ 

+ {a - 6f-^f[-{b + 25) f + g + b\b - 6f-'f] 
-{b- 5f-'f[-ia + 25) f + g + a^a - 5)''~^f]. 

Next, we replace f\'N^b,N-?,,N-2,i^{N-2),^{N-2)\ by (IXSaj) . f\'N^A,N -2,i^{N - 
2),V'(iV-2)| by dAJb]), g\N^A,N -2,i,{N -2),i,{N -2)\ by ([AjED, // by (|ASD, fg by 
()A.3dp and fhhy ()A.3ep . Then we find the remaining is nothing but the l.h.s.of ()A.6p . which is 
zero. □ 

A. 3 Proof in the case of a generic T matrix 

In fact, if tp satisfies ()3.8ap and ()3.8bp . then we can get formulae ()A.2p and we only need a 
modification for the proof of the identity ()A.6p . In the proof we need the following 

Lemma 3. 110 (see also Ulf ) Suppose that 3 is an N x N matrix with column vector set 
{Hj}; 17 is an N X N operator matrix with column vector set {^j} and each entry Qj^s being an 
operator. Then we have 

N N 

^\n,*E\ = Y,m^),*E^l (A.IO) 

i=i i=i 
where for any N 'th- order column vectors Aj and Bj we define 

and 

\Aj *E\ = \Ei, - ■ ■ AjoEj, Ej+i,--- ,En\ 

In the hght of the above Lemma taking ^Ij^g = {E^)^ — 35{E'^)'^ + 35"^ and H = / or / we 
have 

(TV(r)/)/ = (Tr(r)/)/ (A.ll) 
instead of (|A.5p . This equality together with formulae (|A.2p allows us to get the identity 
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